AD-A188  851 


UNCLASSIFIED 


RELATIVE  NOTION  OF  THO  SATELLITES  IN  CLOSE-ELLIPTICAL 
ORBITS(U)  AIR  FORCE  INST  OF  TECH  HRIGHT-PATTERSON  AF8 
OH  SCHOOL  OF  ENGINEERING  J  H  HEINTEL  NOV  87 
AFIT/GA/AA/87D-2  F/G  22/2 


AF I T/BA/ AA/87D-3 


RELATIVE  MOTION  OF  TWO  SATELLITES 


IN  CLOSE-ELLIPTICAL  ORBITS 


THESIS 


James  M  Meintel 
Second  Lieutenant,  USAF 


AFIT/GA/AA/87D-3 


c..;;cTE 


FEB  0  9  1988  r: 


Approved  for  public  release;  distribution  unlimited  ",or1  For 


.a  :  ‘  •  -  M 


vnnnnt 


v 


AF I T/GA/ AA/B7D-3 


RELATIVE  MOTION  OF  TWO  SATELLITES  IN  CLOSE-ELLIPTICAL  ORBITS 


THESIS 


Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 

In  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Master  of  Science  in  Astronautical  Engineering 


James  M.  Meintel,  B.S.E. 
Second  Lieutenant,  USAF 


November  1987 


Approved  for  public  release;  distribution  unlimited 


Pre-face 


The  purpose  of  this  study  was  to  find  the  solution  to 
the  relative  equations  of  motion  for  two  satellites  in 
close-elliptical  orbits.  This  was  done  by  performing  a 
harmonic  analysis  on  the  Floquet  solution,  yielding  an 
expression  for  the  solution  in  terms  of  time  and 
eccentricity. 

The  solution  was  first  verified  for  the  circular  orbit 
case  to  confirm  the  accuracy  of  the  computer  code.  It  was 
then  found  for  the  elliptical  orbit  case  for  small 
eccentrici ties. 

In  finding  the  solution  and  writing  my  thesis,  I 
received  invaluable  support  from  my  faculty  advisor. 

Dr.  William  Wiesel.  Without  his  guidance  this  thesis  would 
have  never  been  completed.  Finally,  I  would  like  to  thank 
my  wife  Laura  and  daughter  Kathryn  for  their  understanding 
during  all  the  extra  hours  that  I  spent  with  this  thesis  at 
their  expense. 


James  li.  Meintel 
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Abstract 


The  need  -for  a  closed  form  solution  for  two  bodies  in 


close  elliptical  orbits  is  identified.  Equations  of  motion 


are  derived  using  Hamilton's  equations.  Floquet  theory,  and 


its  applicability  to  this  problem,  is  discussed.  r  the 


solution  for  circular  orbits  is  derived  in  closed  form  and 


L v^d  in  closed  form  and 

(Tt  -  1  ■  >n  H-  i  O  r  '■  <■'  •  ;  C  /  K-t  "t  >■  '■■■ 


numerically,  using  Floquet  theory  and  harmonic  analysis,  f 


A  <j. ■*'  ,  A/Vy 


Yhe  solution  for  elliptical  orbits  is  found  numerically 


first  for  very  small  eccentricities,  finding  the  dependence 


on  eccentrici ty .  The  solution  is  then  found  for  all  time  as 


a  function  of  eccentricity.  lift  :> 


A *  M  « 


RELATIVE  MOTION  OF  TWO  SATELLITES  IN  CLOSE-ELLIPTICAL  ORBITS 


Introduction 


The  relative  motion  of  satellites  is  a  problem  that  has 
been  widely  studied.  The  results  are  very  important  when 
one  wants  to  rendezvous  or  dock  two  vehicles.  Historically, 
the  first  use  for  a  solution  was  during  the  Gemini  program 
when  the  first  docking  maneuvers  were  perfected  for  use  on 
the  Apollo  missions.  A  complete  solution  to  the  problem  is 
relatively  simple  for  circular  orbits;  however,  elliptical 
orbits  add  a  degree  of  difficulty. 

Numerous  authors  have  published  the  solution  to  the 
problem  of  two  vehicles  in  circular  orbits,  including  Buning 
(1984)  and  Kaplan  (1976).  Buning  derives  the  equations  of 
motion  for  the  elliptical  case,  but  he  limits  solution  to 
the  circular  case.  Investigations  of  the  elliptical  problem 
usually  amount  to  subtracting  the  position  and  velocity 
vectors  of  the  individual  bodies.  This  method  is  limited, 
however,  because  of  a  loss  of  significant  digits  during 
subtraction.  Lancaster  <1970)  formed  a  computational  method 
to  calculate  the  relative  position  and  velocity  for 
elliptical  orbits;  however,  the  results  are  for  single 
points  in  time  and  are  not  a  solution  to  the  equations  of 
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motion.  It  also  does  not  include  out  of  plane  motion. 
Berreen  and  Crisp  <1976)  form  a  solution  for  a  probe  ejected 
into  an  elliptical  orbit  from  a  space  station  in  a  circular 
orbit,  but  does  not  address  the  problem  of  both  bodies  being 
in  elliptical  orbits. 

The  dynamics  for  this  problem  can  be  solved  using  the 
Lagrangian  and  the  Hamiltonian  as  outlined  in  Meirovitch 
(1970).  The  resulting  equations  of  motion  are  linear  and 
periodic,  permitting  a  Floquet  analysis  to  be  done.  The 
system  eigenvectors  and  the  particular  solution  to  the 
equations  of  motion  can  be  expressed  as  functions  of  their 
Fourier  coefficients,  as  outlined  by  Brouwer  and  Clemence 
(1961).  The  results  then  lead  to  a  complete  solution  to  the 
relative  equations  of  motion  expressed  as  an  expansion  of 
the  eccentricity. 
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This  chapter  defines  the  problem,  equations  of  motion, 
and  describes  the  theory  used  to  analyze  the  equations  of 
motion.  The  coordinate  system  is  defined  as  a  rotating 
rectangular  system  with  its  origin  centered  on  body  A.  The 
equations  of  motion  are  derived  using  Lagrange's  and 
Hamilton's  equations.  The  resulting  system  is  linear  and 
time— periodic ,  making  Floquet  theory  applicable  in  solving 
for  the  homogenous  solution  to  the  system.  A  brief 
explanation  of  the  solution  for  autonomous  systems  is  also 
included.  In  addition,  this  chapter  contains  information 
relevant  to  finding  repeated  eigenvectors,  Fourier 
coefficients,  and  a  particular  solution  to  the  equations  of 
mot ion. 


Equations  of  Mot i on 

The  equations  of  motion  are  derived  by  finding  the 
Lagrangian  and  then  forming  the  Hamiltonian  for  the  system 
shown  in  Figure  1.  The  figure  shows  two  coordinate  systems. 
The  first  system  has  its  origin  at  the  center  of  the 
gravi tat i onal  field,  and  the  polar  coordinates  R  and  Q 
describe  the  orbit  of  body  A.  The  vector  q  is  the  position 
of  body  B  with  respect  to  A. 


Figure  1.  Physical  System  and  Coordinate  System 

Both  bodies  are  assumed  to  be  in  two-body,  unperturbed 
orbits.  The  masses  of  each  body  are  also  assumed  to  be 
negligible,  thus  any  gravitational  attraction  between  bodies 
A  and  B  can  be  ignored. 

The  second  coordinate  system  is  a  rotating,  rectangular 
coordinate  system.  The  origin  of  the  system  is  at  the 
center  of  body  A  and  follows  body  A  along  its  two— body 
orbital  path.  The  x-direction  is  along  the  position  vector, 
R,  and  the  y-direction  is  perpendicular  to  the  x-direction 
in  the  orbital  plane.  The  z-direction  completes  the 
right-handed  system.  For  this  system,  as  seen  from  an 
observer  in  body  A,  the  coordinate  directions  are  up. 


■forward  and  left.  The  z-direction  is  the  only  coordinate 


out  of  the  plane  of  the  orbit. 

In  order  to  determine  the  Lagrangian  for  the  second 
particle,  the  position  and  the  square  of  the  velocity  vector 
for  body  B  must  be  found.  The  position  of  body  B  is: 


r=  (R  +  x)  i+yj+zk  (2.1) 

^  ^ 

The  velocity  then,  including  the  rotating  terms,  is: 

r  —  (R  +  x  —  <9  y)  i+(y+0R+0x)  j+zk  (2.2) 


thus , 


r  •  r  =  R2  +  xZ  yZ0Z  +  2Rx-20Ry-20xy  +  yZ 

+  RZ0Z  +  0Z  x 2  +  20xy  +  2R©y  +  2R  02x  +  zZ  (2.3) 

The  gravity  potential  per  unit  mass  of  body  B  for  the 


two— body  problem  is: 


V/m  =  -  p/|r| 


A£ 

R 


-pC(R  +  x)Z  +  yZ+  zZ]~1/Z 


(2.4) 


Expanding  the  denominator  and  dropping  the  terms  of  order 
greater  than  two  yields 


V/m  =  - 

R 


(2.5) 
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With  the  Lagrangian  per  unit  mass 


X  =  T/m  -  V/m 


(2.6) 


and  the  kinetic  energy  per  unit  mass 


T/m  =  -  (r  •  r) 

2  'M  •v 


(2.7) 


the  resulting  Lagrangian  is: 


=  £  r2  +  £  x2  +  £  ©2y2  +  Rx  —  0  R  y  —  0  x  y  +  -  y2 
22  2r  2  7 


+  £  02R2  +  £  02xz  +  0xy  +  0Ry  +  S2R  x  +  £  z2 
2  2  2 


+  fcL  -  JL 


x2 - *£-  y2  -  — tL-  z 2 

R3  2RS  2R3 


(2.8) 


The  Hamiltonian  for  a  system  is  found  from  the 


following  equation: 


H  =  £  p.  q.  -  X 


(2.9) 


were  the  q's  are  the  coordinates  as  defined  before,  and  the 


p's  are  the  conjugate  momenta  which  can  be  found  by 


dX 

p  =  _ _ 


(2. 10) 


The  resulting  momenta  are 


p  =  x  +  R  -  0  y 
x 


(2.11) 


p  =  y  +  0x+0R 
y 


(2. 12) 


p  =  z 
z 


(2. 13) 


After  solving  the  above  equations  for  the  q.  's  and 


substituting  into  Eq  (2.9),  the  resulting  Hamiltonian  is 


. .  1  2  1  2  1  2 
H  =  -  p  +  —  p  +  -p  — 

2  x  2  y  2  2 


+  i  JL 


q2+iJL 

y  2  Ra 


q  -  Rp  +  0  p  q 
z  x  x  y 


-  0  q  p  -©Rp 
x  y  y 


(2. 14) 


To  find  Hamilton’s  equations  of  motion  for  the  p's  and  q's 


the  fallowing  relationships  must  be  used: 


q.  = 


<2.  15) 


Pi  =  " 


(2. 16) 


Substituting  the  Hamiltonian  into  these  two  equations  yields 


the  following  system  of  equations  which  describe  the  motion 


of  body  B  with  respect  to  body  A. 


q  =  p  -  R  +  ©q 
xx  y 


(2. 17) 


q  =  p  -  0q  -  0  R 
y  y  x 


(2.  18) 


qz  =  p^ 


(2.20) 


=  -  H  +  ^  Q  +  g 


—  q  +  0  p 
3  x  y 


(2.21) 


P  = 
y 


—  q  -  0  p 
%3  y  x 


(2.22) 
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U"l 


© 


1 


_  _  p 


3  qz 


(2.23) 


Putting  these  equations  into  matrix  -form  yields  the 
following  set  of  equations. 


•  « 

■ 

“ 

*  •* 

*  •> 

q 

o 

1 

e 

o 

0 

0 

q 

R 

X 

2/u 

X 

_P 

Px 

R3 

0 

o 

6 

0 

0 

Px 

R2 

q 

-  e 

0 

0 

1 

0 

0 

q 

-0R 

y 

. 

”  = 

p 

x 

y 

-  4*  4 

py 

o  - 

e 

R3 

0 

0 

0 

py 

0 

q 

0 

o 

O 

0 

0 

1 

q 

0 

Z 

z 

0 

o 

0 

0 

R3 

0 

Pz 

0 

»  * 

(2.24) 


For  the  case  in  which  body  A  is  in  a  circular  orbit, 
the  above  reduces  to  a  linear,  constant  coefficient  system. 
However,  if  body  A  is  in  an  elliptical  orbit,  R  and  6  are 
time  dependent.  Nevertheless,  both  R  and  B  are  periodic; 
therefore  the  entire  Hami 1 toni an  described  above  is  a 
linear,  time-periodic  system. 


Constant  Coefficient  Systems 


For  body  A  in  a  circular  orbit,  Eq  (2.24)  can  be 
written  in  the  form  x  =  A  x  +  f  where  A  i s  &  constant 

*Nj  ^ 

matrix,  £  is  a  constant  vector,  and  x  is  a  time  dependent 
vector.  As  with  scalar  differential  equations,  x(t>  can  be 
salved  by  finding  the  particular  and  homogenous  solutions 
independently  and  adding  the  solution 


8 


(2.25) 


f) 


3£. 


X  =  X  +  X 
~  ~h  ~p 


The  easiest  of  the  two  solutions  is  to  find  x  .  If  x 

~p  ~p 


is  assumed  to  be  a  constant  vector  then  x  =  0 

~p 


Thus 


Eq  (2.24)  becomes 


0  =  A  x  +  f 

rs#  *vp 


(2.26) 


and  x  is  easily  obtained. 

~p 


On  the  other  hand,  x  is  not  so  easy  to  find.  The 

Mi 


homogenous  differential  equation  is 


-rjr(x  )  =  A  x 
d  t  ~h  ~  ~  h 


(2.27) 


rearranging 


dx 

~h 


=  A  dt 


(2.28) 


integrating  both  sides  yields 


In  (x  )  =  A  t  +  c 


(2.29) 


taking  the  exponent  of  both  sides 


x  =  exp  (At  +  c)  =  c  exp  (At) 


(2.30) 


where  c  is  a  constant  vector  related  to  the  initial 
condi tions. 

Although  this  is  not  a  mathematically  rigorous 
argument,  it  does  show  how  the  solution  is  obtained.  The 


1 


I 


1 


expression  exp(At)  is  similar  to  the  scalar  exponential  in 


that  the  expansion  is 


exp  (At)  =  I  +  At  +  —  (A  t>2  +  -i-  <A  t)3  +  •••  (2.31) 

■*-  ^  ~  2!  ^  3!  ~ 


A*  however  can  be  put  in  the  form 


A  *  F  J  F' 


(2.32) 


where  F  is  the  matrix  of  eigenvectors  of  A  in  column  form. 


and  J  is  a  diagonal  matrix  of  the  eigenvalues  of  A.  J  is 


said  to  be  in  Jordan  normal  form. 


If  it  is  noted  that 


(F  J  F-1) n  =  <F  J  F_1)  (F  J  F-1)  •  •  •  =  F  Jn  F-1  (2.33) 

fV  •>*  »S(  »v  ~ 


then 


2, 


1  ....  a  , 


exp (At )  -  I  +  F  CJtl  F  +  F  C- ,(Jt)‘1  F  +  F  C^(Jt)  1  F 


=  F  exp(Jt)  F" 


(2.34) 


setting  £(t)  =  F  exp(Jt)  F-  ,  the  total  solution  to  the 


differential  equations  of  motion  is 


x(t)  =£(t)  c  +  x 


(2.35) 


Therefore,  finding  the  solution  to  the  equations  of 


motion  for  the  case  where  body  A  is  in  a  circular  orbit  can 


be  done  by  finding  the  eigenvalues  and  eigenvectors  and 


solving  Eq  (2.26). 


The  solution  for  a  constant  coefficient  system  has  been 
well  studied  and  yields  relatively  simple  results.  The 
system  for  this  problem,  however,  is  one  of  time-periodic 
coefficients.  The  solution  for  this  type  of  problem  was 
discovered  by  Floquet  in  the  latter  part  of  the  1800's.  The 
most  common  uses  for  Floquet  theory  is  to  find  the  stability 
of  time-periodic  systems  in  celestial  dynamics.  Few 
studies,  however,  deal  with  finding  a  solution  to  the 
equations  of  motion. 

If  we  start  with  the  system  of  time-periodic 
differential  equations 

x  »  A (t )  x  (2.36) 

where  x  is  the  state  vector,  x  is  the  time  derivative  of  the 
state  vector,  A(t)  is  a  periodic  matrix  with  period  T.  The 
numerical  solution  to  this  problem  is  only  slightly  more 
difficult  than  the  one  for  constant  coefficients  in  A: 
however,  the  process  is  different  so  a  discussion  of  the 
procedure  is  included. 

Since  we  have  a  linear  system,  its  solution  can  be 
described  as 


x  < t )  =  d>  (t,0)  c  ( t  ) 


(2.37) 


where  the  state  transition  matrix,  <£,  has  the  initial 


conditions 

g  (t,0>  =  A ( t )  g  (t,0)  (2.38) 

g  (0,0)  =  £  (2.39) 

where  I  is  the  identity  matrix. 

Floquet  theory  shows  that  g  can  be  -factored  into  two 
matrices  F  and  J,  such  that 

d>  (t,0)  =  F  (t)  exp(Jt)  F-1  (0)  (2.40) 

The  matrix  J  is  a  constant  matrix  most  conveniently  put  in 
Jordan  normal  -form.  The  diagonal  entries  o-f  J  are  the 
Poincare  exponents  which  are  related  to  the  system 
eigenvalues.  The  matrix  F  is  a  time  periodic  matrix  with 
the  same  period,  T,  as  the  original  system. 

For  the  constant  coefficient  system,  F  would  just  be 
the  eigenvectors  of  A.  The  only  difference  between  the 
constant  coefficient  system  and  the  periodic  system  is  that 
F  is  periodic  in  the  latter  case.  Therefore,  solving  the 
Floquet  problem  requires  finding  the  constant  matrix  J  and 
the  periodic  matrix  F  over  a  single  period. 

The  first  step  in  the  Floquet  analysis  is  to  find 
j£(T,0).  This  is  called  the  monodromy  matrix.  The  monodromy 
matrix  is  usually  found  by  numerically  integrating 

g(t,0)  =  A  <t ) - g(t ,0)  (2.41) 
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over  one  period.  Having  <£(T,0)  and  knowing  F(T)  =  F(0) 
results  in 

£(T,0)  =  F  (0)  exp(Jt)  F~‘(0)  (2.42) 

~  ~  >V 

which  can  be  written  as 

exp(Jt)  =  F-1  ( 0 )  g(T,0>  F(0)  (2.43) 

This  shows  that  F(0)  is  the  matrix  of  eigenvectors  for 
g<T,0) . 

Also,  if  X.  are  the  eigenvalues  of  the  matrix  £(T,0>. 
The  Poincar6  exponents,  co.  ,  are  related  to  the  eigenvalues 
by  the  following  relationship 

X  =  exp(u.T)  (2.44) 

t  v 

Where  the  t*>.  are  the  diagonal  elements  of  J,  or 

u>  =  (1/T)  In  (X  )  (2.45) 

i  v 

The  stability  information  for  the  system  is  now 
present.  If  any  of  the  Poincar6  exponents  have  positive 
real  parts,  the  system  is  unstable.  Since  the  stability 
information  is  all  that  is  usually  needed,  this  is  where 
most  Floquet  analyses  stop.  For  this  study,  however,  the 
solution  to  the  equations  of  motion  are  desired;  therefore, 
the  analysis  must  continue. 

In  order  to  find  the  complete  solution,  F(t>  must  be 
found.  Since  it  is  periodic,  F(t)  is  only  needed  over  the 


■first  period.  By  substituting  Eq  (2. 38)  into  Eq  (2.37)  and 


rearranging 

F  ( t )  =  A ( t )  F ( t )  -  F ( t )  J  (2.46) 

where  the  initial  conditions  for  the  equation  is  just  the 

matrix  of  eigenvectors.  So  by  integrating  this  equation 

over  one  period  the  total  solution  can  be  found. 

The  first  problem  that  arises  is  that  F(t)  and  J  can  be 

complex,  making  any  real  analysis  quite  difficult.  They 

can,  however,  be  arranged  such  that  both  are  completely 

real.  F  should  be  arranged  in  column  vectors,  ff,  if  the 

eigenvector  is  completely  real,  then  f  will  simply  be  the 

eigenvector.  On  the  other  hand,  if  there  is  a  pair  of 

complex  eigenvectors  (they  always  appear  in  conjugate 

pairs),  then  the  columns  will  be  f  and  f 

i.  real  l  un 

The  matrix  J  will  no  longer  be  in  Jordan  normal  form, 
but  will  be  in  block  diagonal  form.  Real  go  's  will  remain 
the  diagonal  elements  of  J,  but  the  complex  pairs  of  go  '  s 

~  t 

will  be  appear  in  the  fallowing  diagonal  blocks: 


Re(<o)  Im(w) 

-Im(Go)  Re(c*>) 


(2.47) 


The  matrix  exp(Jt)  is  then  replace  by  the  diagonal 
entries  exp  (cot)  for  the  real  go  ,  and  the  diagonal  blocks 

t  V 


(2. 48) 


0 


exp  ( u> 


CQs(Im(w  )t)  — sin  ( Im(a>  )  t ) ' 

i  t )  l 

v  5in(Ini(u  )t)  cos(Im(co)t> 

L  t 


for  complex  conjugate  pairs. 

In  this  problem,  it  will  also  be  necessary  to  find 


.-i 


F  (t)  over  one  period.  F(t)  is  always  invertable;  however, 
an  easier  method  for  finding  F-1(t)  exists.  If  the  identity 

.-i 


FF  a  I  is  di f f erent i ated  with  respect  to  time  and 
subsituted  from  Eq  (2.46),  the  result  is 


F-1  (t )  =  -  F1  ( t )  A  (t )  +  3  F_1  ( t ) 


(2.49) 


Since  the  problem  that  is  being  dealt  with  is  periodic, 
there  will  be  a  pair  of  repeated  eigenvectors  along  the 
velocity  vector  of  the  two  body  orbit  (Wiesel,  1981).  Since 
the  columns  of  F  must  be  independent,  another  eigenvector 
must  be  found.  If  X  is  the  vector  of  p  and  q.  of  the  two 

~  V  V 

body  orbit,  then  the  repeated  eigenvector,  £^will  be 

dX 


dt 


(2.50) 


and  the  extended  eigenvector  will  be 


dX 


^  = 


dE 


(2.51) 


where  E  is  the  energy  of  the  orbit. 

Now  that  F(0)  and  J  can  be  formed,  Eq  (2.46)  can  now  be 
integrated.  By  using  harmonic  analysis,  as  outlined  in 
Brouwer  and  Clemence,  one  can  find  a  closed  form  expression 
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for  the  elements  of  F.  Since  F  is  periodic  it  can  be 
expressed  as 


F(i  ,  j) 


-  c  +  c  cos  t  +  c  cos  2t  +  • • •  +  -  c  cos  nt 

2  O  1  2  2  n 

+  s  sin  t  +  s  sin  2t  +  • • •  +  s  sin  nt  (2.52) 

12  r. 


where 


2n-l 

c  =  -  y  F(jot)  cos  kja,  k  =  0,1,2,  •  ,n  (2.53) 

k  n  .  Z4 

j=o 

F(jot)  sin  kja,  k  =  0,1,2,-  --.n  -  1  (2.54) 


and 


2  n -1 


-  i  1. 


J  =  1 


a  =  2n/2n  (2.55) 

% 

Therefore,  numerically  finding  n  values  of  F,  spaced  over 
equal  intervals  of  time,  yields  a  complete  solution  for  F. 
Therefore,  a  complete  homogenous  solution  has  been  found. 

Particul ar  Solution 


In  order  to  find  the  complete  solution  to  the  problem  a 
particular  solution,  in  addition  the  homogenous  solution, 
must  be  found.  Eq  (2.24)  is  in  the  form 

x  =  A  x  +  f (t)  (2.56) 


Now,  introducing  the  modal  variables, 


as 


then 


I 


y  =  F-1  < t )  x  +  F  *(t)  x  (2-58) 

substituting  Eq  (2.57)  into  Eq(2.58)  yields 

y  =-F-1  ( t )  A  (t )  x  +  JF-1  ( t )  x  +  F-1  (t )  A  ( t )  x  +  F_1(t)f(t)  (2.59) 

Substituting  Eq  (2.57)  into  Eq  (2.59)  and  simplifying 
results  in 

y  =  J  y  +  F~‘(t)  f(t)  (2.60) 

The  term,  F  *(t)  f (t) ,  is  periodic;  therefore,  it  can  be 
integrated  over  one  period  to  find  the  Fourier  coefficients 
as  was  shown  for  F(t).  The  solution  for  £  can  then  be 
easily  integrated  by  hand  since  it  is  just  a  series  of  sines 
and  cosines. 

Conclusion 


The  geometry  of  the  problem  and  the  equations  have  been 
shown.  A  brief  overview  of  the  methodology  for  solving  the 
equations  of  motion  was  then  presented.  The  following 
chapters  will  utilize  the  procedures  discussed  in  order  to 
find  the  solution  to  the  problem  of  relative  motion  of  two 
satellites  in  neighboring  elliptical  orbits. 


III.  Circular  Sol ution 


Introduction 


This  chapter  will  deal  exclusively  with  the  solution 


for  the  problem  with  body  A  being  in  a  circular  orbit.  It 


will  first  be  solved  as  an  autonomous  system,  and  then  the 


results  will  be  compared  with  the  solution  obtained  by  using 


a  Floquet  analysis  as  outlined  in  the  previous  chapter. 


Autonomous  Sol ution 


The  elements  of  the  A  matrix,  R,  and  0,  for  the 


circular  solution  will  all  be  constant.  Using  canonical 


variables,  R,  /j,  and  0  will  all  be  identically  equal  to  one. 


The  resulting  equations  of  motion  will  be 


0  0  O 


2  0  0 


0  O 


0  0  10  0 


(3.  1 ) 


0-1-1  O  0  0 


0  0  O  0  0 


0000-1  0 


The  first  step  in  finding  the  exact  solution  to  the 


circular  case  is  to  find  the  homogenous  solution  for 


Eq  <3.1).  As  shown  in  the  previous  chapter,  this  is  done  by 


finding  the  eigenvalues  and  eigenvectors  of  the  A  matrix. 


s* 


The  eigenvalues  are  found  by  the  equation 


det (X  I  -  A)  =  0 


(3.2) 


The  resulting  characteri Stic  equation  is 


X4  +  X2  =  0 


(3.3) 


Therefore,  the  eigenvalues  are 


X  =0 

i.2 


(3.4) 


X  =  ■*■  i 

3.4 


(3.5) 


X  =  ♦  i 
3,<S 


(3.6) 


The  eigenvectors  are  then  found  by  the  equation 


(X  I  -  A)  C  =  0 

~  v 


(3.7) 


For  X^  the  resulting  eigenvector  is 


^  =  (  O,  1,  -1,  O,  0,  0) 


(3.8) 


However,  there  is  no  independent  eigenvector  coorespondi ng 
to  X^.  This  vector,  termed  the  generalized  vector,  can  by 
found  using  the  expression  (Reid  1983) 


(A  -  X  I) 


V  £*  -  C* 


(3.9) 


This  yields 


^  =  <  j,  1,  -1,  -  ^  O,  0)' 


<3. 10) 


The  other  -four  eigenvectors  occur  in  complex  conjugate 


pairs.  As  in  the  Floquet  analysis  outlined  in  chapter  2, 


the  eigenvectors  can  be  broken  up  into  their  real  and 


imaginary  parts,  and  J  can  be  put  into  blocks  of  sines  and 


cosi nes. 


The  pairs  of  eigenvectors  are 


^  =  <  |  i,  -  i,  1,  -  i  i,  0,  0) 


^  =  (  O,  0,  0,  0,  1,  i) 


Therefore,  the  fallowing  can  be  written 


2/3  0  -1/2 


1  -1/2 


-1  -1 


0  -1/3  0 


F1  = 


(3. 11) 


(3. 12) 


(3. 13) 


(3. 14) 


exp(Jt)  = 


0  cos  t  sin  t 


0  -sin  t  cos  t 


cos  t  sin  t 


-sin  t  cos  t 


Multiplying  matrices 


—  si  n  t 
2 


-  —  COS  t 

2 


t  +  1  - 


COS  t 


-  —  sin  t 
2 


F  exp(Jt)  = 

*s. 


-1  -  t  -  1 


cos  t 


sin  t 


—  sin  t 
2 


-  —  COS  t 
2 


£<t ,0)  =  F  exp(Jt)  F  = 


2-cos  t 


sin  t 


sin  t 


-2cos  t 


3t-sin  t 


2-cos  t 


1-cos  t 


3t-2sin  t 


-3t+2sin  t  -2+2cos  t  -l+2cos  t  -3t+4sin  t 


-1+CQ5  t 


-sin  t 


-sin  t 


-l+2cos  t 


(3. 15) 


cos  t  sin  t 


-sin  t  cos  t 


(3. 16) 


cos  t  sin  t 


-sin  t  cos  t 


(3. 17) 


So  the  above  matrix  is  the  £  matrix  and  is  the  solution  to 


the  homogenous  set  of  equations. 


Finding  the  particular  solution  requires  finding  the 


solution  to  the  following  matrix  equations 


•  Vt*.  tVi 


;sWW 
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'  0 

1 

1 

o 

0 

0  ■ 

■  0  • 

2 

0 

0 

1 

0 

0 

1 

-1 

0 

o 

1 

o 

0 

X  =  1 

1 

0 

-1 

-1 

o 

0 

0 

~p 

0 

0 

0 

0 

o 

0 

1 

0 

0 

0 

0 

0 

-1 

0 

0 

k.  y 

Resulting  in 


(3. 18) 


x  =  (  O,  1 ,  -1 ,  1  ,  O,  O) 

~p 


(3. 19) 


Putting  the  homogeneous  and  particular  solutions 
together  yields  an  equation  in  the  following  form. 


x(t)  =  0(t,t  )  c  (t  )  +  x 


~p 


(3.20) 


4ft 


Therefore,  the  next  step  in  finding  the  total  solution  is  to 
find  the  constant  vector  c(tQ) .  Evaluating  Eq  (3.20)  at 
time  equal  to  zero  results  in 


x  (0)  =  0(0,0)  c  (O)  +  x 

Aj  ry  A«ip 


Since  £(0,0)  is  the  identity  matrix 


(3.21) 


c  (0)  =  x  (0)  -  x 


(3.22) 


Recalling  Eqs  (2.11),  (2.12)  and  (2.13)  for  body  A  in  a 


circular  orbit  yields  the  following 


X 


X  (t) 


-  y 


+  X  +  1 


(3.23) 


Evaluating  for  time  equal  to  zero  and  substituting  into 
Eq  (3.22)  yields. 


c  (O)  = 


o 

-  y. 


+  1 


y  +  x 
o  o 


-  1 


(3.24) 


i  ■:  i 

Multiplying  out  Eq  (3.20)  with  the  given  results  tor 
g(t,0),  c  (O)  ,  and  x^  directly  yield  the  following  equations 
for  the  position  of  body  B  with  respect  to  A  in  the  original 
coordinate  system  referred  to  in  Figure  1. 


x  (t) 
y  (t) 


z  (t) 


4x  -(3x  +  2  y  )  cos  t  +  x  sin  t 

o  o  'o  o 


(  6  x  +  4  y  )  sin  t  +  2  x  cos  t 
o  o  o 

-  (  6  x  +  3  y  )  t  -  2  x  +  y 
o  '  o  o  o 


z  cos  t  +  z  sin  t 
o  o 


2  y  (3.25) 

o 


(3.26) 

(3.27) 
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The  equations  of  motion  that  have  been  derived  do  not 
directly  give  equations  for  the  velocities,  or  dot  terms. 

Eq  <3.20)  gives  the  equations  for  the  momenta  terms,  and  the 
velocity  equations  can  then  be  found.  The  equations  for  the 
momenta  are 


p  <t)  =  x  (t )  -  y(t) 


=  -  <  3  x  +2y)  sin  t  -  x  cos  t 
o  o  o 


+  <  6  x  +  3y  )  t  +  2  x  -y 
o  o  o  o 


(3.28) 


p  (t)  =  y  +  x  +  1 
y 

=  -  x  sin  t  +  <  3  x  +  2  y  )  cos  t 
o  o  o 

-  2  x  -  y  +1  (3.29) 

o  o 

p  <t)  =  z <t)  =  —  z  sin  t  +  z  cos  t  (3.30) 

I  o  o 

Therefore,  the  solutions  for  the  velocity  components  are 


x  <t) 
y  <t> 
z  <t> 


3  x  +  2  y  +  x  cos  t 
o  o  o 


-  2x  sin  t  +  (6x  +  4y  )  cos  t 

o  o  o 


-z  sin  t  +  z  cos  t 
o  o 


(3.31) 

6x  -  3y  (3.32) 
o  o 

(3.33) 


These  equations  are  the  same  as  derived  in  Kaplan. 
Kaplan,  however,  derived  the  equations  using  Newtonian 
mechanics,  whereas  the  derivation  here  is  using  Hamilton's 
equations. 
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The  exact  solution  derived  above  should  yield  the  same 
results  as  Floquet  theory.  Briefly,  Floquet  theory  consists 


of  the  following  steps 

1.  The  ^  matrix  will  be  integrated  over  one  period. 

2.  The  F  and  J  matrices  will  be  found  at  time  equal  to 
zero  using  the  eigenvectors  and  Poincar6  exponents 
evaluated  from  £  at  one  period. 

3.  F  will  be  integrated  over  one  orbit.  Since  it  is 
periodic,  the  results  are  valid  for  all  time. 

4.  The  fourier  coefficients  for  the  individual  terms 
in  F  are  found. 

5.  The  solution  for  the  g  matrix  for  all  time  is  now 
available  and  the  solution  can  be  found  using  the 
initial  conditions  the  same  as  was  done  for  the 
Hamiltonian  derivation. 

Integrating  £  over  one  period  yields  the  following 
matrix 

'  1 . OOOOOO  0 . 000000 

IB. 849556  1.000000 
-1B.B49556  0.000000 
0.000000  0.000000 
0.000000  0.000000 
0.000000  0.000000 


0.000000 
0. 000000 
1 . 000000- 
0.000000 
0.000000 
o. 000000 


0.000000 
1 8 . 849556 
18.849556 

1 . 000000 
0.000000 
0 .  000000 


0.000000 
0.000000 
0.000000 
0.000000 
1 . 0000O0 
0 . OOOOOO 


0.000000 
0.000000 
0.000000 
0 . OOOOOO 
0 . 000000 
1 . OOOOOO 


Finding  the  eigenvalues  for  this  matrix  is  not 


extremely  difficult,  X  are  all  one.  Finding  the 

l-<5 

appropriate  eigenvectors  is  not  quite  so  easy.  The  first 
problem  that  arises  is  the  that  IMSL  subroutine  eigrf  can 
not  handle  the  problem.  Due  to  the  six  repeated 
eigenvalues,  IMSL  only  finds  three  independent  ei genvectors. 
One  is  a  vector  along  the  velocity  of  body  A  as  expected, 
and  the  other  two  that  it  finds  are  for  the  out  of  plane 
motion.  They  are 

(  0,  1,  -1,  0,  0,  0>T 
(  0,  0,  0,  0,  1,  0)T 
(  0,  0,  O,  0,  0,  1)T 

It  has  already  been  shown  that  one  of  the  eigenvectors 
will  be 


dX 

^  =  df  (3-34) 

where  E  is  the  energy  of  the  orbit  and  X  is  the  position 
vector  of  body  A  with  respect  the  center  of  the  Earth 
evaluated  at  time  equal  to  zero  or  any  integer  multiple  of 
the  period.  The  vector  X,  in  the  coordinates  being  used,  is 


X  =  (  R  ,  0,  0,  v  ,  0,  0) 

~  o  o 

where  v  is  the  linear  velocity  of  body  A. 


(3.35) 

Therefore,  the 


eigenvector  will  be 


O,  o. 


O,  O) 


(3.36) 


dR 


£  =  (  dE  * 


dv 

o 

dE  ’ 


First  R  and  v  must  be  found  as  a  function  of  E. 
o  o 


Starting  with 


E  =  (v  /  2)  -  <  /j  /  R  ) 
o  o 


(3.35) 


and 


v  =  (  /  R  ) 

o  o 


i/2 


(3.36) 


Substituting  Eg  (3.36)  into  Eq  (3.35)  yields 


R  =  - 
o  2E 


(3.37) 


al  so 


v  »  <  -  2  E  ) 
o 


i/2 


(3.38) 


Taking  the  derivative  of  Eqs  (3.37)  and  (3.3B)  with  respect 
to  E  yields  the  eigenvector 


{ '  =  (  ~~  ,  O,  O,  -  (  -  2  E)  1/2,  0,  0)T 


(3.39) 


2E 


Since  canonical  units  are  being  used,  a  =  1  ,  fj  =  1  , 

and  E  =  —1/2  .  Therefore, 


^  =  (  2,  0,  0,  -1,  0,  0 


(3.40) 


Four  of  the  eigenvectors  have  now  been  found;  however, 
two  more  are  still  needed.  Finding  them  by  hand  turns  out 
to  be  fairly  simple  due  to  the  number  of  zeros  present  in  <p. 
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(3.41) 


uni 


WWW 


=  (  o,  1,0,  o,  o,  0)T 

r  =  <  1,  O,  O,  -1,  O,  0)T  (3.42) 

Putting  all  the  eigenvectors  together  gives  the  -following 


F  (O)  » 


O  2 
1  O 
-1  O 
O  -1 
O  O 
O  O 


0  1 
1  O 
0  0 
0  -1 
0  0 
O  0 


O  O 
O  O 
O  O 
O  O 
1  O 
O  1 


al  so 


F-i(0)  = 


o 

1 

o 

-1 

o 

o 


O  -1/3  O 
0  0  1 

110 
0  0-2 

0  0  0 

0  0  0 


O 

o 

O 

o 

1 

o 


o 

o 

o 

o 

o 

1 


(3.43) 


(3.44) 


Since  all  the  eigenvalues  are  equal  to  one,  all  the 
Paincar6  exponents  are  zero.  This  means  that  the  J  matrix 
will  be  zero,  except  an  off  diagonal  one  due  to  the 
generalized  eigenvector  <r  ).  Therefore, 
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0  0  0  0 


0  0  0  0  0  0 


0  0  0  0  0  0 


(3.45) 


0  0  0  0  0  0 


0  0  0  0  0  0 


0  0  0  0  0  0 


1  t  0  0  o  o 


1  0  0  0  0 


exp(Jt)  = 


0  0 


0  0  0 


(3.46) 


0  0  0 


O  O 


o  o  o  o 


0  0  0  0  0 


Now  everything  has  been  found  that  is  needed  to  find 


F(t).  With  F (O)  as  the  initial  conditions,  F  was  integrated 


over  one  orbit.  It  turned  out,  however,  that  the  chosen 


eigenvectors  were  not  periodic.  Instead  of  returning  to 


their  original  values,  the  fallowing  matrix  resulted 


12.6 


F  (T)  = 


-12.6 


(3.47) 


Since  F(0)  *  F(T),  this  is  not  a  periodic  combination 


of  eigenvectors.  However,  due  to  the  fact  that  all  the 


Poincard  exponents  are  zero,  any  linear  combination  of  the 


eigenvectors  is  still  an  eigenvector.  Therefore,  some 
linear  combination  of  the  first  two  eigenvectors  must  be 
found  to  make  the  set  of  eigenvectors  a  periodic  one. 

Evaluating  the  time  derivative  of  F  at  zero  results  in 
the  following 


(3.48) 


What  was  done  to  solve  this  problem  was  to  find  a  set 
of  eigenvectors  that  make  the  time  derivatives  of  the  second 
and  third  row  of  the  second  column  equal  to  zero  at  time 
equal  to  zero. 

The  first  thing  that  was  attempted  was  to  add  a 
multiple  of  to  If 

0  2  0  1  0  0 

1  a  1  0  0  0 

-1  -a  0  0  0  O 

0-1  0-1  0  0 

0  0  0  0  1  0 

0  0  0  0  0  1 


then 


%  % 


(3.51) 


Adding  the  two  matrices  in  this  manner  did  not  have  the 
desired  affect.  The  next  thing  that  was  tried  was  to  just 
have  a  multiple  of  r  .  The  eigenvalue  matrix  at  zero  will 
be 


( 3 . 53 ) 


It  is  easy  to  see  that  if  a  =  3  ,  then  F(0>  will  be  zero 

in  the  desired  positions. 
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The  eigenvector  matrix 


0  2  0 


0  O 


3  0 


0  0  O 


F  (0)  = 


-3  0  0  0  0  0 


(3.54) 


0-1  0-1  O  O 


0  0  0  0 


0  0  0  0  0 


o  -l/; 


o  o 


F_1(0)  = 


(3.55) 


O  -2 


0  O 


was  integrated  around  one  orbit  and  was  periodic. 


The  next  step  was  to  -find  the  Fourier  coefficients  as 


described  in  Chapter  2.  The  results  of  the  harmonic 


analysis  showed  that  each  element  o-f  the  F  matrix  could  then 


be  expressed  as  a  function  of  time  in  the  following  manner 


sin  t 


cos  t 


0  2-cos  t 


sin  t 


F  (t )  = 


O  -2+2cos  t  -2sin  t 


(3.56) 


-1  -sin  t 


-cos  t 


cos  t  sin  t 


-sin  t  cos  t 


Note  that  these  numbers  are  really  accurate  to  ten 


significant  digits.  The  trailing  zeros  were  left  off  to 
make  the  matrices  easier  to  read. 


exp(Jt)  = 


Therefore,  to  find  ^>(t) 


F  (t )  exp(Jt)  = 


2  sin  t  cos  t 
3t  2-cos  t  sin  t 
-3t  -2+2cos  t  -2sin  t 


-1  -sin  t 


-cos  t 


Recalling  that 


0  -l/: 


F  (O)  = 


0  -2 


(3.57) 


cos  t  sin  t I 
-sin  t  cos  tj 
(3.58) 


(3.59) 


0(t,O)  =  F(t)  exp(Jt)  F~  1  (O)  = 

2-2cos  tO  O' 

3t-2sin  tO  O 

3t+4sin  t  0  0 

-l+2cos  tO  O 

O  cos  t  sin  t 

0  -sin  t  cos  t 

(3.60) 

As  expected,  this  £(t,0)  is  identical  (within  ten 
significant  digits)  to  the  one  -for  the  exact  solution.  From 
here  on  the  total  solution  can  be  found  from  Eq  (3.20)  in 
the  same  manner  as  the  exact  solution.  Due  to  redundancy 
this  will  not  be  shown.  However,  the  total  solution  will 
result  in  the  same  solution  because  <g(t,0)  is  the  same  for 
both  approaches. 


2-cos  t  sin  t  sin  t 

3t-sin  t  2-cos  t  1-cos  t 

-3t+2sin  t  -2+2cos  t  - 1 +2cos  t  - 

-1+cos  t  -sin  t  -sin  t 

0  0  0 

0  0  0 


Conclusion 

The  solution  to  the  equations  of  motion  for  the 
circular  case  are  not  new.  The  main  purpose  for  finding 
this  solution  was  to  help  verify  the  computer  programs  that 
were  used.  Since  the  two  solutions  are  identical ,• it  shows 
that  the  numerical  approach  that  was  used  was  correct  for 
the  case  of  a  circular  orbit. 


$ 
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Introduction 


The  procedure  -for  -finding  a  numerical  solution  to  the 
elliptical  problem  is  identical  to  finding  the  numerical 
solution  for  the  circular  case.  The  goal,  however,  is  to 
find  the  relative  motion  solution  as  a  function  of 
eccentricity.  This  will  be  done  by  finding  the  solution  for 
slightly  eccentric  orbits.  For  very  small  eccentricities 

—7 

(ie.  e  =  10  )  only  the  linear  terms  of  e  will  appear. 

After  finding  the  dependence  on  e,  larger  values  for  the 
eccentricities  will  be  used  to  find  the  equations' 

2  3 

dependence  on  e  and  e  . 


Homogenous  Solution  to  the  Equations  of  Motion 


<£<T,0)  is  found  in  the  same  manner  as  for  the  circular 
case.  £  -  A  £  is  numerically  integrated  over  one  period. 
For  all  eccentricities,  g(T,0)  appears  in  the  following  form 


1  0  0  0  0  0 

A  1  0  B  0  0 

C  0  1  D  0  0 

0  0  0  1  0  0 

0  0  0  0  1  0 

0  0  0  0  0  1 

Where  A,B,C,D  are  some  value. 

The  eigenvalues,  as  in  the  circular  case,  are  all  of 
magnitude  equal  to  one;  therefore,  all  of  the  Poincar^ 
exponents  are  zero. 

As  with  the  circular  case,  IMSL  cannot  find  distinct 
eigenvectors  to  the  matrix;  however,  it  is  not  too  difficult 
to  find  the  eigenvectors  by  hand  due  to  the  abundance  of 
zeros  in  the  matrix. 

0  O  O  0  0  O 

A  0  O  B  0  0 

C  O  0  D  O  0 

0  O  0  0  0  0 

0  0  0  0  0  0 

0  O  0  0  0  O 

Four  of  the  eigenvectors  are  easy  to  determine.  They  are 

^  (  0,  1,0,  0,  0,  0)T  (4.3) 


g(T,0)  = 


£  =  (  1,  0,  0,  -A/B,  0,  0)T 

^  (  o,  o,  o,  o,  1,  0)T 

-  <  o,  o,  o,  o.  o.  1)T 


(4.4) 

(4.5) 

(4.6) 
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and,  as  shown  in  Chapter  2,  there  is  a  repeated  eigenvector 
and  its  corresponding  extended  eigenvector.  Recalling 


dX 

s*  -  TF 


<4. 7) 


and 


dX 

^2  =  ~dE~ 


(4.8) 


the  last  two  eigenvectors  can  be  -found.  The  repeated  vector 
will  be 


£*  =  (  O,  p/R*  ,  vo,  0,  0,  O) 


(4.9) 


Where  the  zero  subscript  denotes  time  equal  to  zero.  Using 
canonical  units 


jj/Rz  =  1/  ( 1— e) 2 


(4. lO) 


From  Kaplan  (1976) 


e  = 


(R  v  / 
o  o 


-  1 


(4.11) 


Which  leads  to 


v  =  C  (fj/R  )  (e  +  1)  1 
o  o 


1/2 


(4. 12) 


thus 


v  =  C(e  +  l)/(e  -  1)3 
o 


1/2 


(4. 13) 
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Now  only  the  second  eigenvector  needs  to  be  found 


Recalling  Eq  (3.35) 


E  =  (v  /  2)  -  (  p  /  R  ) 
o  o 


Substituting  for  v  and  solving  for  R  results  in 

o  ^  o 


R  =  <e  -  1)/2E 
o 


thus 


(4. 14) 


(4. 15) 


dR  /dE  = 
o 


1  -  e 
2  E2 


(4. 16) 


However,  for  canonical  units  the  semimajor  axis,  a,  will  be 
equal  to  one,  and  E  =  -1/2  .  Thus, 


B  -V'j 


dR  /dE  =2  <1  -  e) 
o 


Now,  in  order  to  find  dv  /dE,  start  with 

o 


v  =  C  (p/R  )  (e  +  1)  1 
o  o 


1/2 


Using  the  chain  rule 


dv  /dE  =  (6v  /<5R  )  (6R  /dE) 

o  o  o  o 

=  -  (1/4)  C(1  +  e)/R3:1/2  (1  -  e)  /  2E2 

o 


simpl ifying. 


dv  /dE  =  -  C(1  +  e) / ( 1  -  e) ] 
o 


1/2 
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(4. 17) 


(4. 18) 


(4. 19) 


(4.20) 


-.V 


Since  the  eigenvectors  have  all  been  found,  the  matrix 


F (O)  will  be 


2 ( 1-e) 


(1-e) 


F  (0)  = 


Gs]' 


(4.21) 


-m 


-A/B  O  O 


The  next  step  is  to  integrate  F  over  one  period.  This 


was  done  to  the  matrix  above  and,  as  in  the  circular  case, 


it  was  not  periodic.  The  method  for  making  F  periodic  for 


the  circular  solution  was  to  make  the  second  column  of  the 


second  and  third  rows  of  F  equal  to  zero  by  using  a  multiple 


of  the  first  eigenvalue.  However,  there  is  no  constant  that 


will  make  this  happen.  Nevertheless,  if  the  first 


eigenvector  is  multiplied  by  three  as  before,  then 


2(1 -e> 


(1-e) 


F  (0)  = 


GS]' 


(4.22) 


-A/B  0 


w'H  L,-*  _  •  _  *»« 
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F (3,2)  =  O 


(4.23) 


and 


F (2 , 2)  =  -  e/(l  -  e)2 

•v* 


(4.34) 


For  the  case  where  e  =  0  ,  Eq  (4.34)  reduces  to  zero. 

This  method  was  tried  and  resulted  in  a  periodic 
function  for  F(t).  Fifty  values  of  F  were  saved  along  equal 
increments  of  time  during  one  period  and  the  Fourier 
coefficients  were  found  as  discussed  in  Chapter  2. 

Table  4.1  lists  the  eccentrici ties  and  the  calculated 
values  for  the  Fourier  coefficients  of  the  first  column  of 
the  first  row  of  F(t). 


Table  4. 1 

Fourier  coefficients  for  F(l,l) 


e 

S 

s 

s 

1 

•  2 

3 

.0000001 

-.000000300 

* 

* 

.0000010 

-.000003000 

* 

* 

•OOOOIOO 

-.000030000 

* 

* 

.OOOIOOO 

-.000300000 

* 

* 

. 0005000 

-.001500000 

-.000000750 

* 

. 0008000 

-.002399999 

-.000001920 

-.000000002 

.OOIOOOO 

-.002999999 

-.000003000 

-.000000003 

. 0050000 

-.014999859 

-.000074999 

-.000000422 

. 0080000 

-.023999424 

-.000191992 

-.000001728 

.0100000 

-.029998874 

-.000299980 

-.000003375 

. 0200000 

-.059990999 

-.001199680 

-.000026990 

. 0300000 

-.089969627 

-.002698380 

-.000091048 

s  =  -3e  + 

1 

s  =  -3e2 

2 

1 . 125e2 

s  =  -3.375e3 

3 

where  the  stars  indicate  zero  to  nine  significant  digits. 
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V, 

V 

,1 


—7 

For  e  =  10  ,  one  can  see  the  linear  dependence  on 

e.  As  the  eccentricity  increases,  the  square  and  cubic  are 
no  longer  negligible  and  can  be  extracted  -from  the  data. 

The  resulting  equation  tor  F(l,l)  as  a  function  of  both 
eccentricity  and  time  is 

F ( 1 , 1 )  =  <-3e  +  1.125e2)  sin  t  -  3e2  sin  2t  (4.35) 

'v 

—  3.375e3  sin  3t 


The  rest  of  the  data  to  determine  the  coefficients  for 
the  other  elements  of  F(t)  is  in  Appendix  8.  The  resulting 
equations  for  the  rest  of  the  elements  of  F  are  as  follows: 


F(1 ,2) 


F(1 ,3) 


F  ( 1 ,4) 


=  2  +■  e2  +  C-2e  +  (3/4)e33  cos  t  -  e2cos  2t 

-  (3/4)e3  cos  3t 

=  Cl  -  2e  +  (5/8) e2  +  . 75e3!  sin  t 

+  Ce  -  2e2  +  ( 1/3) e3  3  sin  2t 
+  C (9/8) e2  -  2. 25e3 1  sin  3t  +  (4/3)e3  sin  4t 

=  -  e  +  Cl  —  (9/8)eZ3  cos  t  +  Ce  -  (4/3)e3)  cos  2t 

+  (9/8)e2  cos  3t  -  (4/3)  cos  4t 


F  ( 1  ,5)  =  F  ( 1 , 6)  =  0 

F (2, 1 ) ’  =  3  +  1.5e2  +  (6e  +  2.25e3)  cos  t 
+  7.5e2  c os  2t  +  9.75e3  cos  3t 

F (2,2)  =  (-e  +  . 38e3)  sin  t  -  e2  sin  2t  -  (9/8)e3  sin  3t 
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F  < 2 , 3)  =  2 


-  2. 5e  +  e2  -  ,87e3  +  (-1  +  4e  +  3.875e2 


1.5e3)  cos  t  +  (-1.5e  +  5eZ  +  4.16e3)  cos  2t 


<2. 12ez  +  5.8e3)  cos  3t  +  3e3  cos  4t 


F  <2 , 4)  =  Cl  -  (3/8) e2 1  sin  t  +  ( 1 . 5e  -  e3)  sin  2t 


+  2.215e2  sin  3t  -  2.96e3  sin  4t 


F(2,5) 
F (3, 1 ) 


=  F (2,6)  =  0 

•y 


-3  +  1 . 5e2  +  (— 3e  +  1.875e3)  cos  t  -  3e2  cos  2t 


-  3. 375e  cos  3t 


F (3,2) 
F  (3 , 3) 


=  0 


-2  +  . 5e  +  e2  -  .25e3  +  (2  -  2e  -.5e2  +  1.244e3) 


cos  t  +  (1.5e  -  2ez  -  .4e3)  cos  2t 


(1.5eZ  -  2.25e3)  cos  3t 


F (3,4)  = 

"y 


(-2  +  .5e2)  sin  t  +  (-1.5e  +  1.08e3)  sin  2t 


-  1.5e2  sin  3t  -  (5/3>e3  sin  4t 


F (3,5) 
F (4,2) 


=  F (3,6)  =  F(4, 1)  =  O 


-1  +•  .5e2  +  C-e  +  (5/8) e3]  cos  t  -  e2  cos  2t 


+  (9/8) e  cos  3t 


F  ( 4 , 3 )  = 


C-l  +  (7/8)e23  sin  t  +  C-e  +  (7/6)e3]  sin  2t 


-  (9/8) e2  sin  3t  +  (4.3)e3  sin  4t 


F(4,4)  =  C-l  -  (3/8)e2)  cos  t  +  C-e  -  (l/6)e3]  cos  2t 


-  (9/8)e2  cos  3t  -  (4/3) e3  cos  4t 


it.  «».  a*. 


•M 


?] 


^.5>w 


F (5,5)  =  -1 . 5e  -  1  - 5e  -  1 . 5e  +  Cl  +  e  +  <5/B)e  +  <5/B)e  3 

cos  t  +  ( .  5e+  .5eZ  +  . 168e3)  cos  2t 
+  C(3/B)e2  +  (3/8)e33  cos  3t  +  (l/3)e3  cos  4t 
F  (5,6)  =  Cl  -  e  -  (l/8)e2  (l/8)e33  sin  t  +  ( .  5e  -  .  5e2 

-.165e3)  sin  2t  +  C(3/8)e2  -  (3/8)e33  sin  3t 


+  (l/3)e3  sin  4t 


F  (6, 1 )  =  F (6,2)  =  F (6,3)  =  F (6,4)  =  0 


F  (6,5)  = 


F (6,6) 

"y 


=  C— 1  -  e  -  (5/8) e2  -  (5/8)e33  sin  t 

+  C-e  -  e2  - ( 1/3) e3  3  sin  2t 
+  (~1.125e2  -  1.125e3)  sin  3t  -  1.35e3  sin  4t 
=  1  -  e  -  (l/8)e2  +  <l/8)e3  +  Ce  -  e2  -  (l/3)e33cos  t 
+  C(9/8)e2  +  (9/8)e33  cos  3t  +  (4/3)e3  cos  4t 

(4.36) 


Since  F(t)  has  been  found,  finding  <£(t,0)  is  now 
relatively  simple.  As  in  the  circular  solution 


0  0 


exp ( Jt )  = 


(4.37) 


Finding  the  solution  for  F-1(0>  is  tedious,  but 
relatively  simple.  It  will  be  a  matrix  in  the  form 


1-V  v ->Vv 


» 

I 
♦ 

’^•VC 


<1-1 f )  /c 


-g/c 


F  <t>  = 


(4.38) 


where 


=  3/(1  -  e) 


=  -  3  t  <  1  +  e)  /  <  1  -  e)  3 
=  2  <1  -  e) 

=  -  C  (1  +  e)  /  <  1  -  e)  31/2 


=  -  A/B  (from  Eq  (4.1)) 


=  -d/ <ec  -  d) 


=  c/<ec  -  d) 


(4.39) 


The  total  solution  for  £<t,0)  can  now  be  expressed  in  terms 
of  eccentricity  cubed  by  the  following  matrix  multiplication 


£(t,0)  =  F  < t )  exp  < Jt )  F  (0) 

^4  r^» 


(4.40) 


Combining  Eqs  (4.36),  (4.37),  and  (4.38)  the  homogenous 

solution  for  x  =  A(t)  x  is  now  determined  as  an  expansion 

»S>  *\i 

of  eccentricity  in  Eq  (4.40). 


Parti cul ar  solution 

The  complete  solution  for  the  relative  motion  is  almost 
complete.  All  that  is  needed  is  the  particular  solution. 
Recalling  Eq  (2.59) 

y  =  J  y  +  F -1  (t )  f(t)  (4.41) 

Since  F  4(t)  f (t)  is  periodic  it  was  integrated  and  the 
Fourier  coefficients  were  found  as  done  for  F(t).  The 

solution  for  ^  is  then  easy  to  integrate  by  hand  and  results 
in  the  following: 

y  ( 1 )  =  K  +  K  +  t  — 1  +  (2/3) e  -  (l/3)e2  +  -5e3]  t 
12 

+  C (4/3)  -  (S/3)  +  (7/6>e2  -*•  (l/3)e3!  sin  t 
+  Ce  -  (ll/6)e2  +  .4125e3D  sin  2t 
+  C (9. 43) e2  — 1 . 5e3 3  sin  3t  +  e3  sin  4t 

*(2)  =  Kz 

^(3)  =  K  +  < -e  -  . 75e2 )  .t  +  (-2  -  1.75e2)  sin  t 

+  C- (3/2) e  -  (2/3) e3  3  sin  2t  -  (17/12)e2sin  3t 

-  1.48e3  sin  4t 

^(4)  =  -  (2  -  2. 75e2)  cos  t  -  C(3/2)e  -  2.4125e3)  cos  2t 

—  (17/12)eZcos  3t  -  1 . 48e  3cos  4t 

^(5)  =  £(6)  =  O  (4.42) 

Since  the  particular  solution  can  be  any  function  that 
satisfies  the  differential  equation,  any  initial  value  for  ^ 


will  be  valid;  however,  ^(0)  =  <-•  leads  to  easier  expressions. 
This  can  be  done  by  settinq  K  =  O  and 

1-3 

K  ■  2  +  (3/2) e  +  (4/3) e2  +  0.9325e3  (4.43) 

4 

Now  >£(t)  has  been  found,  the  particular  solution  is 

x  =  F ( t )  y  ( t )  (4.44) 

■ — p  ~ 

where  F(t)  is  expressed  in  Eq  (4.36)  and  ^(t)  in  Eq(4.42) 


>1  ete 


sol ut i on 


Recalling  Eq  (2.35), 

x(t)  =  #<t,0)  c  +  x  (4.45) 

~ 

however,  since  ^(O)  =  0 

c  =  x  (4.46) 

~  ~o 

and  Eq  (4.45)  can  be  written 

x(t)  =  d>(t,0)  x  +  x  (4.47) 

where  g(t,0)  is  expressed  in  Eq  (4.40),  x^  in  Eq  (3.23),  and 
x  in  Eq  (4.44) 


Cone 1 us i on 


Using  the  expressions  derived  in  this  chapter,  the 
solution  to  the  equations  o f  motion  can  be  expressed  as  a 
■function  o-f  eccentricity  and  time.  The  results  reduce  to 
the  circular  solution,  but  the  actual  equations  for  the 
elliptical  solution  are  rather  complicated  to  write  out  in 


an  expanded  form. 


V 


Conclusion 


An  expression  for  the  solution  to  the  relative  motion 
problem  has  been  found.  Although  this  solution  is  limited 
to  small  eccentricities  (the  error  is  approximately  e3 ) ,  it 
has  wide  applications.  It  could  be  applied  to  traffic 
management  around  the  space  station.  Shuttle  activities,  or 
any  other  application  that  utilizes  orbits  of  small 
eccentri ci ties. 

These  equations  of  motion  are  an  improvement  over 
previous  methods.  They  are  not  limited  to  body  A  being  in  a 
circular  orbit  as  in  some  previous  studies'.  The  results  are 
valid  even  at  very  close  distances,  which  is  a  weakness  of 
the  methods  based  on  subtracting  the  two  position  vectors. 

This  analysis  also  has  the  advantage  of  being  able  to 
calculate  the  necessary  changes  in  velocity  for  rendezvous. 
If  we  recall  the  basic  solution  to  the  equation  of  motion 

x  (t)  =  pit  ,0)  x  +  x  (5.1) 

The  values  of  ^(t,0),  x^ ,  and  x^  have  previously  been  found. 
For  the  rendezvous  problem,  the  initial  and  final  times  are 
known.  Also,  the  initial  position  is  known,  and  the  final 
position  is  simply  zero  for  the  three  different  directions. 
The  only  values  that  are  not  known  are  the  initial  and 
final  velocities.  This  results  in  a  system  of  six 
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equations  and  six  unknowns.  These  are  not  difficult  to 


solve  since  all  the  other  values  are  constants. 

One  difficulty  in  using  this  solution  occurs  when  the 
desired  initial  time  is  not  at  perigee,  but  at  some  time,  t  . 
This  problem  can  be  overcome  by  using  the  fact  that 

g(t,0)  =  ^(t,to>  g(to,0>  (5.2) 

Finding  the  required  solution,  therefore,  requires  solving 
Eq  (5.2). 

There  are  two  limitations  to  this  solution.  The  first 
is  that  the  relative  distances  between  the  two  bodies  must 
be  small.  The  second  limitation  is  that  the  solution  is  for 
small  eccentr i ci ti es.  Equations  for  greater  eccentricities 
can  be  found  using  the  same  approach  as  this  thesis.  This 
could  be  done  with  the  same  computer  code,  finding  the 
dependence  on  higher  orders  of  eccentricity. 

Utilization  of  these  equations  yield  themselves  best  to 
computer  analysis  since  they  are  relatively  lengthy 
expressions.  The  equations  being  somewhat  unwieldy  are  only 
a  minor  inconvenience.  One  must  remember,  however,  that  the 
equations  are  in  canonical  units. 
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Introduction 


This  appendix  is  meant  to  describe  the  code  that  was 
used  for  this  thesis.  Hopefully,  this  will  make  it  easier 
to  understand  the  numerical  processes  used. 

Computer  Code 

The  main  programs  were  relatively  simple.  Their 
purpose  was  to  initialize  some  of  the  parameters,  mainly 
eccentricity,  and  call  the  numerical  integrator.  There  were 
actually  three  main  programs.  The  first  integrated  the  £ 
matrix  over  one  orbit  and  its  output  was  just  £(T,0).  The 
second  integrated  the  F  matrix  over  one  period  and  its 
output  was  the  values  of  the  F  matrix  for  fifty  evenly 
spaced  time  intervals  over  one  period.  The  last  integrated 
F  1  over  one  period  and  the  output  was  fifty  evenly  spaced 
values  of  F  1  f (t),  where  f(t>  is  the  forcing  function. 

The  numeric  integrator  that  was  used  was  Haming.  It  is 
a  fourth  order  predi ctor-corrector  capable  of  integrating 
systems  of  first  order  differential  equations.  Using  it 
entails  forming  a  subroutine  called  "rhs,"  where  rhs 


calculates  the  right  hand  side  of  the  equations  of  motion. 


The  position  o-f  body  A  was  necessary  in  the  subroutine 
rhs,  so  it  was  calculated  using  a  Newton-Raphson  method  to 


■find  Kepler's  equation. 

The  fourier  coefficients  were  found  using  the  approach 
defined  in  chapter  2. 
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Appendix  B:  Data  From  Harmonic  Anal vses 


Table  B. 1 

Fourier  Coefficients  for  F(l,2> 


e 

C 

o 

c 

i 

c 

2 

c 

3 

. OOOOOO 1 

2.00000000 

-.000000200 

* 

* 

.0000010 

2.00000000 

-.000002000 

* 

* 

.0000100 

2 . 00000000 

-.000020000 

* 

* 

.0001000 

2.00000001 

-.000200000 

-.000000010 

* 

. 0005000 

2.00000025 

-.001000000 

-.000000250 

* 

. 0008000 

2.00000064 

-.001600000 

- . 000000640 

* 

.OOIOOOO 

2.00000100 

-.002000000 

- . 00000 1 000 

* 

. 0050000 

2.00002500 

-.009999906 

-.000025000 

-.000000094 

. 0080000 

2.00006400 

-.015999616 

-.000063997 

-.000000384 

.0100000 

2.00010000 

-.019999249 

-.000099994 

-.000000750 

. 0200000 

2.00040000 

-.039994000 

-.000399894 

- . 000005998 

. 0300000 

2.00090000 

-.059979751 

-.000899460 

-.000020233 

c  =  2  +  e2 
o 

c  =  -2e  +  (3/4) e9 

i 

2 

C  —  —<B 

2 

c  =  — (3/4) e9 

3 

Table  B.2 

Fourier  Coefficients  for  F(l,3) 


.0000001 
.OOOOOIO 
.0000100 
.0001000 
. 0005000 
. 0008000 
. OO 10000 
. 0050000 
. 0080000 
.0100000 
. 0200000 
. 0300000 


. 999999800 
. 999998000 
. 999980000 
. 999800006 
.999000016 
. 998400400 
. 998000626 
.990015719 
. 984040383 
. 980063246 
. 960255944 
. 940582466 


s 

=  1  -  2e 

l 

s 

=  e  -  2e 

2 

S 

=  (9/8)  e 

3 

S 

=  (4/3)  e 

000000010 

000000100 

000001000 

000099980 

000499500 

000798720 

000998000 

004950042 

007872176 

009800347 

019202878 

028210067 


. _ 3 


00000001 1 
000000281 
,  000000719 
000001 123 
,  000027844 
,  000070848 
,000110251 
,000432018 
,000951858 


000000001 

000000165 

000000672 

000001307 

000010239 

000033836 


e  +  (5/8)  e  +  <3/4)e 


>> 


Table  B.4 

Fourier  Coefficient  for  F(2,l) 


e 

c 

o 

c 

i 

c 

2 

c 

a 

.0000001 
.0000010 
.0000100 
.0001000 
. 0005000 
. 0008000 
.0010000 
. 0050000 
. 0080000 
.OlOOOOO 
. 0200000 
. 0300000 

3. 00000000 
3.00000000 
3.00000000 
3.00000002 
3.00000038 
3.00000096 
3.00000150 
3.00003750 
3.00009600 
3.00015001 
3.00060018 
3.001 3509 1 

. 00000060 
. 00000600 
. 00006000 
. 00060000 
. 00300000 
. 00480000 
. 00600000 
.03000028 
.04800115 
. 06000225 
. 12001800 
. 18006080 

* 

* 

* 

. 000000075 
.000001875 
. 000004800 
. 000007500 
.000187500 
. 000480004 
.000750010 
.003000160 
.00675081 1 

# 

* 

* 

* 

. 00000000 1 
. 000000005 
.000000010 
. 000001219 
. 000004992 
. 000009750 
. 000077996 
. 000263222 

c  =  3  + 
o 

L.5eZ 

c  =  6e  + 

1 

2.25  e3 

c  =  7.5e2 

2 

c  =  9. 75e3 
a 

Table  B.5 

Fourier  Coefficients  for  F<2T2> 


e 

s 

s 

s 

1 

2 

a 

.0000001 

-.000000100 

* 

# 

.0000010 

-.000001000 

* 

* 

.0000100 

-.000010000 

* 

* 

.OOOIOOO 

-.000100000 

-.000000010 

* 

. 0005000 

-.000500000 

-.000000250 

* 

. 0008000 

-.000800000 

-.000000640 

* 

.0010000 

-.001000000 

-.000001000 

- . OOOOOOOO 1 

. 0050000 

-.004999953 

-.000025000 

- . OOOOOO 1 4 1 

. 0080000 

-.007999808 

-.000063997 

-.000000576 

.0100000 

-.009999620 

-.000099992 

-.000001124 

. 0200000 

-.01 9996997 

-.000399894 

-.000008997 

. 0300000 

-.029989877 

-.000899461 

-.000030350 

s  =  -e  + 

1 

2 

s  =  -e 

2 

.  38e3 

s  =  -1. 125e3 
a 
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Fourier  Coefficient 


Table  B.7 

Fourier  Coefficient  for  F<2,4> 


e 

S 

s 

s 

s 

O 

l 

2 

a 

.0000001 

1 . 000000000 

.000000150 

* 

* 

.0000010 

1 . 000000000 

.000001500 

* 

* 

.0000100 

1 . ooooooooo 

. 000050000 

* 

* 

. 0001000 

. 999999996 

. 000150000 

. 000000021 

* 

. 0005000 

. 999999906 

. 000750000 

. 00000053 1 

* 

. 0008000 

. 999999760 

.001200000 

.000001360 

. 000000002 

.0010000 

. 999999625 

.001499999 

.000002125 

. 000000003 

. 0050000 

. 999990625 

.007499885 

.000053129 

. 000000370 

. 0080000 

. 999976000 

.011999531 

. 000135993 

.000001515 

. 0 1 00000 

. 999962500 

.014999084 

.000212482 

. 000002958 

. 0200000 

. 999849996 

.029992666 

.000849711 

. 000023656 

. 0300000 

. 999662488 

.044975252 

.00191 1039 

. 000079808 

s  =  1-  .375e2 

1 

s  =  1 . 5e 

2 

a 

-  e 

s  =  2. 125e2 

3 

s  =  2.96e3 

4 

Table  B.8 

Fourier  Coefficient  for  F(3,l) 


e 

c 

o 

c 

l 

c 

2 

c 

3 

.0000001 
.0000010 
.0000100 
.0001000 
. 0005000 
. 0008000 

. 0010000 
. 0050000 
. 0080000 
.0100000 
. 0200000 
. 0300000 

-3.00000000 
-3. 00000000 
-3.00000000 

-2.99999999 
-2.99999963 
-2.99999904 
-2.99999850 
-2.99996250 
-2.99990400 
-2.99985000 
-2 . 99939994 
-2.99864970 

-.000000300 

-.000003000 

-.000030000 

-.000300000 

-.001500000 

-.002400000 

-.002999998 

-.014999766 

-.023999040 

-.029998124 

-.059984996 

-.089949371 

* 

* 

* 

-.000000030 
-.000000750 
-.000001920 
-. 000003000 
-.000074998 
-.000191990 
-.000299975 
-.001199600 
-.002697975 

* 

* 

* 

* 

* 

- . 000000002 
-.000000003 
-.000000422 
-.000001728 
- . 000003375 
-.000026989 
-. 000091037 

■N 


CM  N  M  O  O' 
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O  CM  CO  M3  O'  CO 
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000000 
000000 
000000 
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Table  B. lO 

Fourier  Coefficient  for  F(3,4) 


e 

S 

1 

5 

2 

5 

a 

s 

4 

.0000001 

-2.00000000 

-.000000150 

* 

* 

.0000010 

-2.00000000 

- . 00000 1 500 

* 

* 

.0000100 

-2.00000000 

-.000015000 

* 

* 

.0001000 

-2.00000000 

-.000150000 

-.000000015 

* 

. 0005000 

-1.99999988 

-.000750000 

-.000000375 

* 

. 0008000 

- 1 . 99999968 

-.001199999 

-.000000960 

* 

.0010000 

-1.99999950 

-.001499999 

- . 00000 1 500 

- . 000000002 

. 0050000 

-1.99998750 

-.007499865 

-.000037499 

- . 000000208 

. 0080000 

-1.99996800 

-.011999445 

-.000095993 

-.000000853 

.0100000 

-1.99995000 

-.014998917 

-.000149984 

-.000001667 

. 0200000 

-1.99980006 

-.029991334 

-.000599737 

-.000013326 

. 0300000 

-1.99955033 

-.044970761 

-.001348672 

-.000044943 

s  =  -2  + 

1 

.  5e2 

s  a  -1.5e  +•  1.08e3 

a 

s  =  — (5/3) e3 

4 

Table  B. 1 1 

Fourier  Coefficient  for  F(4,2) 


e 

c 

o 

c 

l 

c 

2 

c 

a 

.0000001 

-1.00000000 

-.000000100 

* 

♦ 

.0000010 

-1.00000000 

-.oooooiooo 

* 

* 

.0000100 

-1.00000000 

-.000010000 

* 

* 

.0001000 

-0.99999999 

-. 000100000 

-.000000010 

* 

. 0005000 

-0.99999988 

-.000500000 

-.000000250 

* 

. 0008000 

-0.99999968 

-.000800000 

-.000000640 

* 

.0010000 

-0.99999950 

-.000999999 

-.000001000 

-.000000001 

. 0050000 

-0.99998750 

-.004999922 

-.000024999 

-.000000141 

. 0080000 

-0.99996800 

-.007999680 

-.000063997 

-.000000576 

. OlOOOOO 

-0.99995000 

-.009999376 

-.000099991 

-.000001125 

. 0200000 

-0.99979998 

-.019995000 

-.000399866 

- . 000008996 

. 0300000 

-0.99954990 

-.029983123 

-.000899325 

-.000030346 

c  =  -1  +  .5ez 
o 

c  =  -e  +  . 625e3 

1 

2 

c  =  -e 

2 

c  =  1 . 125e3 
a 
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Table  B. 12 

Fourier  Coefficient  for  F(4,3) 


.0000001 
.0000010 
-OOOOIOO 
.0001000 
. 0005000 
. 0008000 
.0010000 
. 0050000 
. 0080000 
.0100000 
. 0200000 
. 0300000 


-1 . 00000000 
-1.00000000 
- 1 . 00000000 
-0.99999999 
-0.99999978 
-0.99999944 
-0.99999913 
-0.99997813 
-0. 99994400 
-0.99991250 
-0.99965001 
-0.99921257 


,000000100 
.000001000 
,  OOOO 10000 
.000100000 
,  000500000 
,  000799999 
,  000999999 
.  004999854 
,  007999403 
.  009998833 
.019990668 
,029968508 


* 

* 

* 


,00000001 1 
,  00000028 1 
,  000000720 
,000001125 
,000028124 
,  000071993 
.000112484 
.000449741 
,  001011 191 


* 

* 

* 

* 

* 

* 


OOOOOOOOl 
000000167 
000000683 
00000 1 333 
000010659 
000035945 


s  = 


s  = 


-1  +  . 
-e  +  ( 
-(9/8) 
(4/3) 


875e 


7/6) e* 
2 


Table  B. 13 

Fourier  Coefficient  for  F(4,4) 


e 

c 

1 

c 

2 

c 

3 

c 

A 

.0000001 

-1.00000000 

-. 000000 10O 

* 

•* 

.0000010 

-1.00000000 

- . 00000 1 000 

* 

* 

.0000100 

-1 . 00000000 

- . OOOO 1 OOOO 

* 

* 

.0001000 

-1.00000000 

-.000100000 

-.OOOOOOOl 1 

* 

. 0005000 

-1.00000009 

-. 000500000 

-.000000281 

* 

. 0008000 

-1.00000024 

-.000800000 

- . 000000720 

* 

.001 OOOO 

-1.00000037 

-.001000000 

-.000001125 

-.OOOOOOOOl 

. 0050000 

-1.00000938 

-.005000020 

-.000028125 

-.000000167 

. 0080000 

-1.00002400 

-. 008000085 

-. 000072000 

-.000000683 

.0100000 

-1.00003750 

-.010000167 

-.000112499 

-.000001333 

. 0200000 

-1.00015005 

-.020001334 

-.000449989 

-.000010665 

. 0300000 

-1.00033776 

-.030004506 

-.001012443 

-.000035990 

c  =  -1  -  (3/8) e2 

1 

c  =  -e  -  (l/6)e3 

2 

c  =  - (9/8) eZ 

3 

c  =  - (4/3) e9 

A 
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Table  B. 15 

Fourier  Coefficient  for  F(5,6) 


e 

s 

1 

s 

2 

s 

a 

s 

4 

.0000001 

. 999999900 

. 000000050 

# 

* 

.0000010 

. 999999000 

. 000000500 

* 

* 

.0000100 

. 999990000 

. 000005000 

* 

* 

.0001000 

. 999899999 

. 000049995 

. 000000004 

* 

. 0005000 

. 999499969 

. 000249875 

. 000000094 

* 

. oooaooo 

.999199920 

. 000399680 

. 000000240 

* 

.0010000 

. 998999875 

. 000499500 

. 000000375 

* 

. 0050000 

. 994996891 

.002487479 

. 000009328 

. 00000004 1 

. ooaoooo 

.991992064 

.003967915 

. 000023807 

.000000169 

.0100000 

. 989987625 

.004949835 

.000037123 

. 000000330 

. 0200000 

.979951001 

.009798693 

.000146967 

. 000002612 

. 0300000 

. 969890879 

.014545635 

. 000327209 

. 000008724 

s  =  1  -  e  -  <l/8)e2  H 

1 

2 

s  =»  .  5e  ~  .  5e  -.165c 

2 

s  =  (3/8)  e2  -  (3/8)  eS 
9 

s  =  (l/3)e3 

4 

►  (l/8)e3 
a 

» 

Table  B. 16 

Fourier  Coefficient  for  F(6,5) 


B 

s 

1 

5 

2 

5 

a 

s 

4 

.OOOOOOl 

-1.00000010 

- . 000000 1 00 

* 

* 

.OOOOOIO 

-1.00000100 

-.000001000 

* 

* 

.OOOOIOO 

-1.00001000 

-.000010000 

* 

* 

.0001000 

-1.00010001 

-.000100010 

-.OOOOOOOl 1 

* 

. 0005000 

-1.00050016 

-.000500250 

-.000000281 

* 

. 0008000 

-1.00080040 

-. 000800640 

-.000000721 

* 

.0010000 

-1.00100063 

-.001001000 

-.000001126 

-.000000001 

. 0050000 

-1.00501570 

-.005025042 

-.000028266 

-.000000167 

. 0080000 

-1.00804032 

-.008064172 

-.000072576 

-.000000688 

.0100000 

-1.01006313 

-.010100337 

-.000113626 

-.000001347 

. 0200000 

-1.02025511 

-.020402721 

-.00045901 1 

-.000010879 

. 0300000 

-1.03057992 

-.030909281 

-.001042934 

-.000037073 

s  =  -1  -  e  -  (5/8) e2 

1 

s  =  -e  -  e2  -  (l/3)es 

2 

s  =  - (9/8) e2  -  (9/8)* 
a 

s  =  - (4/3) e3 

4 

-  (5/8) e3 

a 

a 

'-■.v.v. 


WWUJW1WWWBUWJUUUUWWUUNVUP 


uiviuniP  YWTWWXTFTX 


wwwwwnwimtT 
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Table  B. 17 

Fourier  Coefficient  for  F<6,6) 


e 

c 

1 

c 

2 

c 

3 

c 

4 

.0000001 
.0000010 
.0000100 
.0001000 
. 0005000 

.0010000 

. 0050000 
. 0080000 
.0100000 
. 0200000 
. 0300000 

. 999999900 
. 999999000 
. 999990000 
. 999899999 
. 999499969 
. 998999875 
. 994996890 
.991992064 
. 989987625 
.979951001 
. 969890879 

. 000000100 
.000001000 
.000010000 
. 000099990 
.000499750 
. 000899000 
.004974959 
.007935831 
. 009899670 
.019597387 
.029091271 

* 

* 

* 

.00000001 1 
.000000281 
.000001124 
. 000027984 
.000071421 
. 0001 11369 
. 000440901 
.000981628 

* 

* 

* 

* 

* 

. 00000000 1 

. OOOOOO 1 66 
. 000000677 
. 000001320 
.000010450 
. 000034895 

c  =  1  -  e  -  <l/8)eZ  +  <l/8)e3 

1 

c  =  e  -  e2  -  <l/3)ea 

2 

c  =  <9/8>e2  +  <9/8)e3 

9 

c  =  (4/3)ea 

4 
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O  O  O  (0  in  N  O' 

ooooo^n 
*  *000000-1 
o  o  o  o  o  o  o 
o  o  o  o  o  o  o 
o  o  o  o  o  o  o 

_ r  r  i*  r  r  i*  r 

NNOOOMnfflW 

t-owininNOM> 

OOSNNOOOH 

o-(Nin>OMinoin 

OOOOOOCMOCM 

0000-0<fMD 

ooooooo-«(o 
o  o  o  o  o  o  o  o  o 
ooooooooo 


o  o  o  o  o  ^ 
o  o  o  o  o  o 

10  O  O  O  O  O 

o  io  o  o  o  o 
0  0(0000 
oooMint 

O  O  O  O  -«  CM 

o  o  o  o  o  o 


HNNMOO 
O  *0  O  10  IN  M 

O  ^  O  10  >0  o 

O  O  M  — 1  O  >0 
OONO-M 

o  o  cm  o  o  o 
io  in  o  o  o  o 


oooooooooooo 

r  r  i*  r  r  r  r  r  r  r  r  r 


COON^N 

oooo<r~ 
ooooo- 
o  o  o  o  o  o 

o  o  o  o  o  o 

o  o  o  o  o  o 

o  o  o  o  o  o 

o  o  o  o  o  o 

a  a  a  a  a  a 

CM  CM  CM  CM  CM  CM 

I  I  I  I  I  I 


o  o  o  o 
o  o  o  o 
o  o  o  o 
in  co  o  o 
o  o  -  in 
o  o  o  o 
o  o  o  o 


in  in  -  -«  cm  O' 

rv  N  O  O  CM  o 

(0  CM  in  O  MJ 

O  ^  N  o  N 

o  o  — •  s  in 

o  o  o  o  o  -« 

o  o  o  o  o  o 

o  o  o  o  o  o 

■  a  a  ■  a  a 

CM  CM  CM  CM  CM  CM 
I  I  I  I  1  1 


«  nJ.nt 


l_*  *4  -  J 


,  iV 


.‘Al/iL'n  '»r  i»<  i 


,1  »-4  , 


Table  B.20 

Fourier  Coefficient  for  y(4) 
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2.00000000 

. 000000300 

* 
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2.00000000 
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* 
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2.00000000 

. 000030000 

* 

* 
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. 000000042 

* 
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1.99999931 
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.000001062 
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. 002399998 
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0010000 

1 . 99999725 

.002999995 
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. 000000006 
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1.99993125 

.014999396 

.000106245 

. 000000740 

0080000 

1 . 99982400 
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.000271968 

. 000003029 
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.000005915 

0200000 
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0300000 

1 . 99752558 
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.003818637 

.000159449 
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